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SUMMARY  To  opt  to  use  geometrical  methods  to  design  actual 
nonlinear  systems,  one  must,  first  of  all,  resolve  the  problems  of  the 
complexity  of  transformations,  the  great  amount  of  calculations,  and 
other  similar  problems.  Because  of  this,  simplifying  the  design 
process  is  extremely  important.  In  this  article,  based  on 
transformational  ideas,  we  have  obtained  two  methods  and  used  them  in 
the  design  of  the  U.S.  F-8  Crusader  fighter  plane's  vertical  control 
■system-  The  simulation  tc-ults  clearly  demonstrate  that  the  new 
control  laws  are  obviously  superior  to  the  original  nonlinear 
optimization  control  laws.  Moreover,  it  is  possible  to  guarantee  that 
aircraft  can  make  high  angle  attack  flights. 


Key  Words  transformation,  nonlinear  system,  feedback  design 

I.  INTRODUCTION 

Do  to  the  excellence  of  linear  system  design  tools.,  at  the 
present  time,  relevent  nonlinear  systems  design  theory  depends  Co- 
help  in  every  case  on  the  design  methods  of  linear  systems.  In  the 
theory  of  differential  geometry,  tc ansf ormation  methods  are  based  on 
this  idea.  The  key  to  this  is  nothing  other  than--in  looking  for  a 
type  of  transform  —  taking  the  original  system  transform  t  ;>  b?  a  1  i  i  ?  \  ■ 
system  or  quasi!  inea-  system  in  a  different  space.  After  dj/io^  that , 
it  is  possible,  in  the  new  space,  to  opt  for  the  ao  of  linei-  s/s  ten 
iasign  tools  to  oar-y  out  the  design.  Because  almost  all 
transformations  om  '•>  -  -j;ogni/.ad  as  differential  homeomorph  isms ,  and 
differential  homeomorphism  guarantees  that  the  st-ictural 
..'hi- icter  istics  of  systems  before  and  after  transformations  is 
i nvar iibUj ,  iin  this  way, going  again  through  reverse  transformations, 
one  obtains  the  o-ijinal  system  control  laws  which  can  cause  the 
original  system  pole  points  not  to  follow  the  operating  points  in 
their  changes.  T’v-*  difficulty  with  this  method  is  in  getting  th3 
transformational  relationships.  Actual  or  practical  transformation 
f  >'  ns  general  have  a  base  transformat  ion  and  a  feedback 

i.r  .  ns  for  mat i on .  Tr  ansf  or  m  i 1  i  on  methods  can,  thus,  be  of  various  types 
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and  kinds. 

There  are  those  systems  whose  nonlinear  natures  are  due  to  an 
inappropriate  selection  of  base.  Naturally,  it  is  possible  to  opt  for 
the  use  of  overall  transformations ,  taking  them  and  changing  them  into 
standard  forms  of  integral  patterns  that  can  be  controlled.  Hunt  and 
others  in  Reference  (1)  ,  on  the  basis  of  solutions  to  linear  partial 
differential  equations,  gave  out  sufficient  conditions  for  overall 
transformations  as  well  as  general  methods  of  constructing  them. 
Polynumeric  actual  or  real  systems  are  only  capable  of  partial 
transformation  into  linear  systems.  Su  R.,  in  Reference  (2),  proved 
or  confirmed  the  sufficient  conditions  for  the  existence  of  partial 
Inferential  homeomorphisms .  As  far  as  the  size  of  the  partial  domain 
for  this  origin  point  is  concerned,  there  still  is  no  conclusion. 

Fur  thermo'* e ,  the  general  run  of  systems  are  only  capable  of 
transforming  into  linear  systems  within  a  tangent  space.  In  this 
regard,  there  is  Baumann  along  with  others  who  brought  up  extended 
linearization  and  Reboulet  and  others  who  brought  up 

B166 

pseudolinearization  .  The  two  methods  are  basically  in  agreement 
with  each  other.  Even  if  on i  gets  the  closed  ring  character istic 
values  set  up  when  the  closed  ring  operations  are  at  different 
equilibrium  points  to  be  invariable,  because  thn  is  so,  on-  avoids 
point  to  point  linearization. 

F:'om  Reference  (5),  it  is  possible  to  see  that  general  forms  of 
high  order  system  transformation  relationships  are  always  extremely 
difficult  to  obtain.  There  are  only  block  triangular  form  systems 
that  are  an  exception.  Because  of  this,  designers,  when  designing 
iutonatic  aircraft  control  systems,  always  take  flight  equations  and 
do  simplification  processing  in  order  to  obtain  block  triangular 
forms.  However,  this  necessarily  carries  with  it  a  certain  error. 

What  follows  deals  with  investigations  into  the  F-3  fighter  plane's 
vertical  control  system  and  how  transformations  were  selected  fo*  use 
to  arrive  at  control  rules. 
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THE  SET  UP  OF  VERTICAL  MOVEMENT  EQUATIONS  FOR  AIRCRAFT 


II  . 


Modern  generations  of  high  performance  aircraft  normally  fly 
within  large  envelopes  or  ranges,  and  nonlinear  factors  severly  impact 
the  dynamic  responses  of  systems.  A  common  example  is  high  angle  of 
attack  flight  in  fighter  planes.  Their  lift  coefficients  are 
certainly  not  capable  of  using  linear  angle  of  attack  functions  for 
their  precise  expression.  At  high  angular  velocities  of  roll  and  turn 
nonlinear  inertial  combinations  lead  to  loss  of  stability.  Tn  order 
to  improve  dynamic  response,  it  is  necessary,  on  the  basis  of 
nonlinear  models,  to  carry  forward  the  design  work. 

In  Reference  6,  one  takes  lift  coefficients  and  uses  third  degree 
angle  of  attack  functions  to  simulate  actual  lift  curves.  Cosine 
functions  take  the  initial  two  terms  of  a  Taylor  expansion,  at  the 
same  time  taking  into  consi deration  all  the  many  nonlinear  factors, 
under  conditions  in  which  aircraft  have  no  acceleration,  the  Mach 
number  is  0.85,  the  altitude  of  flight  is  9000m,  and  the  condition  is 
level  flight.  On 2  goes  th-ough  a  trimming  down  process  and  obtains 
'  ■  a  jrs  cions  for  the  vertical  movement  of  the  F-8  fighter  plane  (let  the 
-’ugle  of  at  feck  a  -x  ^  let  the  pitch  angle  0  =x  2>  let  the  pitch 
velocity  6  =  .  and,  let  th-  rjdder  deviation  angle  be  6  )  . 


f'  — 0.877*, +  0.47xf  +  3.84 6*f— 0.019*1  —*’*,—  0.088*,*,  +  *,  'j 

i, 

- 

*, 

“*3, 

.  —  4.208*,  — 0.47*,’  — 3. 564*?  — 0.396*,  J 

'  —  0.215  ' 

. 

+ 

0 

t. 

>-20.967. 

In  Reft-rnce  (6),  in  consideration  of  the  linear  optimization 
control  law  (  <5  =  g  )  : 


-0.053*, +  0.5*,  +  0.521*, 
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When  the  initial  perturbation  value  for  the  angle  of  attack 
rtt.>0.43rad(2 1.5')  the  system  loses  stability.  Because  of  this, 
one  finds  the  third  order  form  for  nonlinear  optimization  control 


/*  =  /'»+  0.04*}  —0.048*i#»  +  0.374*f  —  0.312*!#i  ^  ^ 

(2)  This  equation  comes  from  solving  Hamilton-Jacobi  equations. 
The  equation  in  question  can  only  guarantee  $  (V  is  the  Lyapunov 
function)  in  its  negative  specification.  It  does  not  guarantee  the 
positive  specification  of  V.  Because  this  is  true,  when  the  initial 
perturbation  value  a0 >0 . 48r ad ( 27 . ^  )  ,  (1)  is  still  unstable.  Below, 
we  opt  for  the  use  of  a  new  method  of  transformation  to  solve  for  even 
more  satisfactory  results. 


III.  TWO  TYPES  OF  TRANSFORMATION  METHODS 


1 .  Method  1 

Consider  the  system 

i(f> -/<*«» + 


(3) 


In  this,  analyzing  /,  giR"  'a  the  domain  or  region  adjacent  to  the 
origin,  /(0)  =  0,  It  can  be  seen  that  f  and  g  make  the  vector  field  e  *° 
including  the  origin  open  up  on  R  "*  .  Our  objective  is  to  search  for 
the  linear  independent  tr ans for mat  ion  T(x): 

R’-*Rn  that  will  take  system  (3)  and  change  it  into  the  type  of 
integral  form  of  quasicontrollable  standard  form. 


Tt  > 

fi 

T , 

+ 

0 

J 

i 

• 

• 

• 

t.  > 

1 

B 1  67 
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In  this,  Tn+1  should  be  the  arbitrary  function  Tj/...,Tn, 
is  easy  to  see  that 


It 


6T 

6x 


9T  f.9T 


Taking  the  equation  above  and  comparing  it  with  (4),  it  is  possible  to  get 


~^-g=  (dT,,g)  >=  0  »'“1, 2, •••»*— 1 

Therefore,  one  has 

<dT„g)  =  ld<dT„f),g) 

=  -~-(«d'/,g)  =  0 

(dTt-ifg)  =  (  — tg)—Q  (5) 

idTm,gi  =  =  l 


In  this  (ad1  f  ,g ) =<<ad1 f ,g>  ,g> 

Let  the  controllable  array  /  =  [g  ,( ad1 f  ,g  ),...,( adn  f,g)]. 

Then,  equation  (5)  can  be  generalized  to  be 

-g£-/  =  CO,0,~,0 ,1] 

0% 


If  /  is  nonsingular;  select  /*'  's  final  line.  Assume  it  is 
q(x)  (q(x)  is  a  lxn  dimensional  vector  quantity).  Then,  it  is 
possible  to  obtain 


dTi/dx=q(x) 


(6) 


If  q(x)  is  a  gradient  field,  then  (6)  is  a  soluble  equation.  If 
this  is  not  the  case,  then,  normally,  it  is  possible  to  place  a 
function  c-U"  causing  ?(*)«=*(*)•  q(x)  to  be  a  gradient 

field.  Then,  from  {8Tl/dx=q(x)  it  is  possible  to  solve  for 

( x )  .  At  this  time,  system  (3)  transforms  to  become  like  the  form 
below 


( 

ti 
•  • 

T. 


+ 


C 

0 


« 


n  ■ 


in.,  this  method  .and  using  it  on  system  (1),  it  is  first 
/  lu  solve  the  controllable  ar-nay  / 


j  =  (_g,(.ad'  f  ,g),  ( adlf,g )3 


In  this 


( ad'f,g)  =  -df/dx-g 


31  68 
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After  that,  it  is  still  necessary  to  solve  for  d~x  .  Also,  solving 
equation  (6),  the  amount  of  calculations  can  be  imagined.  In  a  practical 
study  of  system(l):  the  influence  of  i>  on  the  system  is  only  -0.019 x^_ 
of  the  a  signal.  Due  to  th:  fact  that  0  is  a  slowly  transf owning 
modality,  its  connection  with  and  influence  on  the  quickly  changing 
modalities  of  «  and  i  are  very  small.  Because  this  is  th»  case,  we 
would  be  well  served  if  we  first  carried  out  designing  of  the  interior 
ring  composed  of  the  quickly  changing  modalities 


rii  a  r  — 0.877*1  +  0.47*1  +  3.846*!  -*!*i-0.088*i*i  +  xt 
*•  i.J  k- 4. 208*. -0.396*,-0. 47*! -3. 564*! 


-0.216  a 

k 

-20.967J 


r  /• 

/. 


vr0'215  >. 

J  1 — 20.967J 


(7) 


System  (7)  's  Jacob  array  or  matrix  is 


F  = 


3/ 

dx 


I"  —  0.877  +  0.94*,  +  1  1 .538*!  —2*,*,  — 0.088*3 
k-4. 208-0. 94*,  —  10.7*J 


Let 


(ad'  f  ,g)  =  % 


fa .(*)  q 

^<*i(*)  J 


+  0.088*i-*! 

-0.396 


In  this 


So  W  i  ng 


a,(x)  =  0.21 5Fn  +  20. 967F„,  «,<*)  =0.215fsi  +  20.967F« 

d  ' s  final  line,  one  obtains 


Taking 


q(x)  20.967a,(*)-0.215a,(*) 

«?  (x  '-2J.  967-a  (x)-0.215a 


■C20. 967-0. 2151 

2(x) 
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it  ’ 3  possible  to  obtain  T j- 20 . 967x ^-0 . 21 5x ^ 


Because 


7\=  =~-f  =  2Q.  967/.-0. 215/, 

f,=  (20.967F,,-0.215F„)i;  +  (20.967F„-0.215F«)i, 


From  T  2>  it  is  possible  to  obtain 

T,~  (20.967F,|-0.215Fii)/i  +  (20.967Ft,  -0.215F«)/» 


Since  this  is  true,  system  (7) 's  status  feedback  rule  is 


4,=(_61r1-6JrI-r,)/ c-»(*)D 

=  (b,T  i  +  biTt  +  Fi)/^(*) 

In  this,  b1  a,ni  b  0  a  r  e  close  :  •'  .  n  j  typ  ?  mu  1 1  i  te-‘  n  equation 
#h  o'  a  :!.•.*•'  I  -  ■;  ;  ;>- f  f  i;  i-nts  ,  that  is,  when  #i  and  *i 

closed  rir.c:  pole  points,  one  has  ) (s  —  e2)  =sI  +  b  s  +b-. 


2.  Method  2 

Consider  the  second  o-nier  affine  nonlinear  system 


Uj  l/,d)J  lf,(D 


] 


w 


It  is  always  capable  of  to'  ansfo»  r>:  i  rig  i  nfo  a  block  tiidnguiar 
Recais-,  in  the  most  yen-”- a  1  situation,  when  y  (  x  )  •  j0  (  \  )9t0  . 

•  >  p  ia  i  i  on  (9)  in  be  written  to  be 


(3) 


r-  e 


(9) 


one . 


8 


l  gl«  H  /,*,  Hg,giy 
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Let  ii  =g  ^-g  2>  Then,  equation  (10)  can  be  written  to  become  a 
block  triangle  form.  As  far  as  its  overall  or  global  transformation  is 
roncerned,  it  is  possible  to  let 

f  1=2/ T9  1^2 


The  n  , 


T'~)  Sid* ,  -j  gldx, 


Because 


2  T  ,  f  [g  2-  f  2g  L 


.  h~  r  fore 


T  2~  P  (  x  ) '  x  1+q(x).x 


In  t n  i  s 


,  ,-x  d/<  f  ,  /■  -f.AlL 

d*i  '  1  3*i  '  9*i 

_v  _  9/l,  +  /■  liL^ALLg  _  /■ -®£i 

?(,)"  5*:  g,  +  /l  3*2  9*;*’  /,3iJ 


If  one  writes  it  to  become  a  controllable  standard  form  of 
expression,  then  one  has 


Tt(x)  =  />(*)  •/,  +  y(x) .  /, 
*(*)  =  - (/>(*)£■ +  ?(*)•*,) 
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In  the  same  way,  it  is  possible  to  ootain  the  status  feedback  rule  or 
law 


u=(.-blTl-b,Tt-Tt)/(-o(x)) 


(ID 


-f  one  takes  system  (7)  and  substitutes  into  the  various 
equations  discussed  above,  it  is  possible  to  see  that  equation  (11) 
and  equation  (8)  are  of  the  same  type. 

IV.  VERTICAL  FEEDBACK  COMTROL  LAWS  AND  SIMULATION  RESULTS 

As  far  as  system  (1)  is  concerned,  due  to  the  fact  that,  after 

the  placement  of  the  pole  points  of  small  sealed  rings,  there  is  no 

relationship  with  the  operating  points,  as  a  result,  it  is  necessary 

to  stabilize  the  integral  signal  6  o£  *  •  It  is  only  necessary, 

after  the  control  law  of  small  sealed  rings  to  attach  a  linear 

f eedb.3c^  K  t?  .  As  a  result  of  this,  one  obtains  the 

(unclear) 

vertical  control  law 


8,  =  {bi'f  i  +  biTi  +  Ti)  I  o  (x)  +  Kt  0 


(12) 


In  this  Tj,  T  2>  T^,  and  «(*)  have  definitions  which  can  be 
seen  in  the  sections  above. 


If  one  takes  the  interior  ring  pole  points  .nd  positions  them  at 
(-2,0)  and  (-4,0),  and,  again,  on  the  basis  of  the  quality  of  actual 
and  simulated  responses,  selects  K(unciear)  =0#*'  and'  below,  if  one 
takas  control  .law  (2)  and  control  law  (12)  and  makes  use  of  them  to 
make  a  complete  comparison  of  the  responses  after  the  original  system, 
the  results  of  simulations  are  as  shown  below: 


1.  When  the  original  perturbation  a0  =0 . 25r*d  ( 14 . 3°  )  ,  the 
responses  of  the  two  angles  of  attack  a,  are  closely  in  line  with 
each  other,  however,  due  to  the  fact  that  the  initial  rudder  offset 
angles  were  not  the  same,  control  law  (12)  was  slightly  better  (see 
Fig.  1  )  . 


]  0 


2.  When  a,  =0 . 47rad  (  27°  )  ,  the  response  of  control  iaw  (12) 
is  clearly  superior  to  control  law  (2)  (see  Fig. 2). 


Wi)WG»d> 


Fig .  1 

A  Comparison  of  Responses  When  ®*  =0.25  rad  (1) 
Law  (2)  (2)  Solid  line  control  law(12) 


F  i  g  .  2 

A  Comparison  of  Responses  When  ®»  =0.47rad  (1) 

Control  Law  (2)  (2)  Solid  line  Control  Law  (12) 

3.  After  >0.47rad,  control  law  (2)  begins 

to  diverge.  However,  straiqht  up  to  the  time  when 
a  _  0.6rad(3l  .4' )  the  response  of  control  law  (12)  is 
( see  F  ig . 3 )  . 


i  I 


Dotted  line  control 


Dotted  line 

to  cause  systems 
still  good 
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Fig. 3  Responses  of  Control  Law  (12)  When  a,  =0.6rad 


Fig.  4  A  CoTp--*ison  of  Rudder  Offset  Angle  Response  Curves  (Curve  1 
represents  the  rudder  offset  for  control  law  (12)  when  a,  =0.47rad. 
Curve  2  represents  the  rudder  offset  for  control  law  (2)  when 
«»  =  o.47ra<'  Curve  3  represents  the  rudder  offset  for  control  law 

(12)  when  «•  =0.6rad) 


V.  SEVERAL  POINTS  IN  CONCLUSION 

This  article,  based  on  the  tr ansformation  thinking  in 
differential  gentry,  designed  for  the  vertical  control  systems  of 
the  F-8  fighter  plane  nonlinear  computer  control  laes  or  rules.  Th 
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results  of  simulations  clearly  show  that  these  control  laws  are 
capable  of  guaranteeing  that  the  aircraft  can  fly  at  high  angles  of 
attack.  The  dynamic  qualities  are  much  better  than  those  of  the 
original  nonlinear  optimization  controls.  A  practical  analysis  of  the 
rudder  offset  responses  in  Fig. 4  is  capable  of  showing  one  that,  when 
angles  of  attack  are  relatively  large,  due  to  6  C^/  <0  (C^  is 

the  lift  coefficient),  there  are  problems  with  the  aircraft  fighting 
the  controls.  However  ,  control  law  (12)  still  actualizes  feedback  to 
defeat  this  before  it  happens,  and,  as  a  result,  it  avoids  this 
fighting  of  the  controls. 

Speaking  in  terms  of  aircraft  in  general,  the  vertical  direction 
is  due  to  the  seperate  characteristics  of  0  and  o.  ,  as  well  as  4  . 
Because  this  is  true,  the  two  methods  that  have  been  discussed  above 
both  are  feasible.  As  far  as  high  order  nonlinear  systems  are 
concerned,  if  it  is  possible  to  set  up  laws  to  take  them  and  transform 
them  into  block  triangular  forms,  then,  both  of  the  two  methods  are 
made  simpler  and  easier. 
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